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ABSTRACT 

We report the most accurate X-ray position of the X-ray source in the giant globular 
O ' cluster Gl in M31 by using the Chandra X-ray Observatory, Hubble Space Telescope 

(HST), and Canada-France-Hawaii Telescope (CFHT). Gl is clearly detected with 
^ . Chandra and by cross-registering with HST and CFHT images, we derive a Icr error 

^ ' radius of 0'.'15, significantly smaller than the previous measurement by XMM-Newton. 

We conclude that the X-ray emission of Gl is likely to come from within the core radius 
P\J ■ of the cluster. We have considered a number of possibilities for the origin of the X-ray 

' emission but can rule all but two scenarios out: it could be due to either accretion 

■"nI" , onto a central intermediate- mass black hole (IMBH), or an ordinary low-mass X-ray 

'nI" ■ binary (LMXB). Based on the X-ray luminosity and the Bondi accretion rate, an 

IMBH accreting from the cluster gas seems unlikely and we suggest that the X-rays 
are due to accretion from a companion. Alternatively, the probability that a 1.5M0 
cluster LMXB lies within the 95 per cent X-ray error circle is about 0.7. Therefore we 
cannot rule out a single LMXB as the origin of the X-ray emission. While we cannot 
. distinguish between different models with current observations, future high-resolution 

and high-sensitivity radio imaging observations will reveal whether there is an IMBH 
at the centre of Gl. 
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1 INTRODUCTION 

Intermediate-mass black holes (IMBHs) have been a sub- 
ject of debate for a long time. If they exist, IMBHs rep- 
resent the long sought after link between stellar-mass and 
super-massive black holes. Until recently, we only have had 
indirect evidence for IMBHs via X-ray observations. For in- 
stance, ultraluminous {Lx > W^^ ergs^^) or hyperluminous 
{Lx > lO^'^ergs"^) X-ray sources have been the best can- 
didates (e.g. Wolter et al. 2006; Farrell et al. 2009) based on 
their X-ray luminosities which, assuming isotropic emission, 
are far in excess of the Eddington limit for stellar mass black 
holes. Furthermore, X-ray spectroscopy and X-ray timing 
also provide some support for IMBHs (e.g. Miller et al. 2004; 
Kong & Di Stefano 2005; Strohmayer & Mushotzky 2009). 



* E-mail: akong@phys.nthu.edu.tw 

t Kenda Foundation Golden Jade Fellow 



However, we still require dynamical evidence in order to con- 
firm solidly that IMBHs exist. 

The only dynamical mass measurement claimed for an 
IMBH candidate is the globular cluster Gl in M31. Gl is 
the most luminous star cluster in the Local Group, and also 
one of the most massive at (4 — 7) x 10^ Mq (Ma et al. 
2009; Barmby et al. 2007). Based on Keck and Hubble Space 
Telescope (HST) observations, it has been claimed that Gl 
hosts a ~ 2 X 10* Mq object at the core (Gebhardt et al. 
2002,2005). However, the suggestion of an IMBH is contro- 
versial and has been challenged by Baumgardt et al. (2003). 
Recently, X-ray emission near the core of Gl has been dis- 
covered based on XMM-Newton observations (Trudolyubov 
& Priedhorsky 2004; Pooley & Rappaport 2006; Kong 2007) 
and it is suggested that the X-rays come from Bondi accre- 
tion from cluster gas onto a central IMBH. However, it is also 
possible that the X-ray emission is due to an ordinary low- 
mass X-ray binary (LMXB), or a collection of faint LMXBs 
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near the core. By refining the relative astrometry between 
XMM-Newton and HST data, Kong (2007) concluded that 
we cannot distinguish between these two scenarios. As we 
will show in this Letter, if there exists an IMBH in Gl, ac- 
cretion likely comes from a companion star (see, e.g. Patruno 
et a. 2006 for a scenario with a giant). 

Alternatively, radio observations may be able to provide 
additional information about the nature of the X-ray source 
in Gl. Ulvestad et al. (2007) employed the Very Large Array 
(VLA) to obtain a low-resolution (~ 3" beam size at 8.4 
GHz) image of Gl and a radio source was detected near the 
XMM-Newton source, about an arcsecond from the cluster 
core. The radio/X-ray flux ratio of Gl (~ 5 x 10~^) is at 
least a few hundred times higher than that expected for a 
LMXB near the cluster centre (~ 5 x 10~*; see Fender & 
Kuulkers 2001), but is consistent with the expected value 
for accretion onto a 2 x W^Mq IMBH (Merloni et al. 2003). 
Furthermore, the radio/X-ray flux ratio is much lower than 
that of supernova remnants (~ 10~^), but it is consistent 
with a pulsar wind nebula (Ulvestad et al. 2007). 

However, it is still a question whether the X-ray and 
radio emission come from the same source because the res- 
olution of both observations do not have sufficient accuracy 
to determine the precise position. While the VLA is still in 
its short baseline configuration and MERLIN is being up- 
graded, the first step is to use Chandra to obtain an accurate 
position for the X-ray source. 

In this Letter, we localized the X-ray emission of Gl by 
performing precise relative astrometry using Chandra, HST, 
and Canada- France-Hawaii Telescope (CFHT). In Section 
2, we describe our X-ray and optical observations and data 
analysis. We present the localization of the X-ray source in 
Section 3. We finally discuss the nature of the X-ray source 
in Section 4. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Chandra 

We observed Gl with the Chandra X-ray Observatory on 
2008 September 30 for a total exposure time of 35 ks (Ob- 
sID 9525). The observation was taken using the Advanced 
CCD Imaging Spectrometer array (ACIS-S) with the tele- 
scope aim point at Gl. Data were telemetered in the very 
faint mode and were collected with a frame transfer time 
of 3.2 s. We used CIAO version 4.fl ACIS ExtraclQ, and 
XSPEC version 12.fjf| packages to perform data reduction 
and analysis. We reprocessed the raw data to make use of 
the very faint mode. In order to reduce the background, only 
events with photon energies in the range of 0.3-7.0 keV were 
included in our analysis. We also inspected the background 
count rates from the SI chip and no flaring event was found 
in the data set. In this paper, we only consider the S3 chip 
of ACIS-S. 

Discrete sources in the Chandra images were found 
with wavdetect (Freeman et al. 2002) together with expo- 
sure maps. We performed source detection on the 0.3-7 keV 
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image. We set the detection threshold to be 10~ , corre- 
sponding to less than one false detection due to statistical 
fluctuations in the background. We performed source detec- 
tion using sequences of wavelet scales that increased by a 
factor of y/2 from scales 1 to 16. A total of 28 X-ray sources 
were detected. The X-ray source in Gl was clearly detected 
with 126 counts. 

We extracted the energy spectrum from a 2" circular 
region centred on Gl. For the background, we selected a 
source- free region with a radius of 15". Response matrices 
were generated by CIAO. We then fltted the background- 
subtracted spectrum with an absorbed power-law model. In 
order to employ statistics, the spectrum was grouped 
into at least 10 counts per spectral bin. Since the statis- 
tics in the spectrum are poor, and therefore it can be fit 
equally well with numerous spectral models (e.g. power-law, 
bremsstrahlung, and thermal plasma) , although a blackbody 
model can be ruled out. However, the absorbed power-law 
model is adopted as it is typical for a globular cluster X-ray 
source. The power-law model provides the best fit {x^/dof = 
7.56/9) to the data with iV^ = 4.71^;? x 10^° cm"^ and 
r = l-Stol (90 per cent confidence level). The 0.3-7 keV 
unabsorbed luminosity is (2.31q'5) ^ 10^^ erg s~^ assuming a 
distance of 780 kpc (Macri et al. 2001). The X-ray spectrum 
and luminosity are very typical for a globular cluster X-ray 
source in M31 (Kong et al. 2002a; Di Stefano et al. 2002). 

The X-ray luminosity as well as the spectrum of Gl are 
consistent with previous XMM-Newton observations (Pooley 
& Rappaport 2006) indicating it is a persistent source. We 
also calculated the Kolmogorov-Smirnov statistic to examine 
whether Gl is a variable source during our Chandra obser- 
vation; there is over 25 per cent probability that Gl is a con- 
stant source. We therefore cannot reject the null hypothesis 
that Gl is a constant source. We also applied the Gregory- 
Loredo variability algorithm (Gregory & Loredo 1992) in 
CIAO to examine our Gl data and the variability index is 
0, implying no variability. 

2.2 Optical Observations 

Gl was observed with the //ST Wide Field Planetary Cam- 
era 2 (WFPC2) on 1995 October 2 with a total integration 
time of 30 minutes in the F814W filter. We downloaded 
the F814W image from the Hubble Legacy Archiv43 for 
which cosmic-ray free, science-quality images are dithered, 
co-added, and corrected for astrometry. By comparing to 
the 2MASS catalogue (Skrutskie et al. 2006), the astrome- 
try of the image is better than 0'.'2. Although Gl was also 
observed with the Advanced Camera for Surveys (ACS) in 
High Resolution Channel (HRC) mode (see Kong 2007), we 
do not use the ACS/HRC data because they have a much 
smaller field-of-view (29" x 26") which is not useful for re- 
fining the relative astrometry with ground-based telescope 
and Chandra data. 

Since the fleld-of-view of WFPC2 is small compared to 
Chandra, it is difficult to improve the relative astrometry be- 
tween HST and Chandra. We therefore obtained a wide-field 
optical image centred on Gl with the MegaPrime/MegaCam 
at the CFHT on 2008 September 5. The MegaCam has an 
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Figure 1. i?5T WFPC2 F814W images of Gl. Tiie red circle is the 95 per cent error circle (0''36 in radius) of the Chandra source and 
its centre is marked by a plus sign. The white dash circle at the centre is the core radius (0'.'2; Ma et al. 1997; Barmby et al. 1997) of 
Gl, while the white circle is the 95 per cent error circle (1'.'47 in radius) of the VLA source. 



array of 36 CCDs, giving a total of 1 degree by 1 degree 
field-of-view. We obtained a series of i'-band images under 
a seeing condition of C/.'? with a total exposure time of 721 
seconds. The images were processed by Terapi^ll which pro- 
vides coadded and astrometrically calibrated images for data 
analysis. By matching over 3000 stars in the field with the 
2MASS catalog, the astrometry of the CFHT image is ac- 
curate to about 0'.'3. 



3 LOCALIZATION OF THE X-RAY SOURCE 
IN Gl 

To localize the precise position of the X-ray source in Gl, 
we have to improve the relative astrometry of the X-ray and 
optical images. To achieve this, we need to match the Chan- 
dra, HST, and CFHT images into a common reference frame 
(see, e.g. Lu et al. 2009). We first compared the Chandra 
image with the CFHT image. By comparing the Chandra 
source list with the CFHT image, we found 5 likely matches 
(in addition to Gl) that are optically bright and isolated. 
They are likely to be foreground stars or background galax- 
ies. Based on these 5 matches, we corrected the astrometry 
of the Chandra image using the IRAE0 task ccmap. The re- 
sulting registration errors are 0'.'062 in R.A. and O'.'ll? in 
declination. 

We next transformed the astrometry of WFPC2 to the 
CFHT image. We found 55 stars in both images that ap- 
peared stellar and unblended. We obtained an astrometric 
solution giving residuals of 0'.'037 in R.A. and 0'.'036 in decl. 
After registering all the images to the CFHT image, we lo- 
cated the optical centre of Gl in the WFPC2 images and 
the position of the X-ray source. We determined the cen- 
troid of Gl in the WFPC2 image (R.A. =00h32m46. 532s, 
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decl.=-)-39d34m40.58s with la error of 0'.'002) by fitting the 
elliptical isophotes to the globular cluster using the IRAF 
task ellipse. 

For the X-ray position, we used ACIS Extract to 
derive the mean position of events within the 90 per 
cent point-spread-function in the 0.3-7 keV band, yielding 
R.A. =00h32m46. 532s, decl.=-|-39d34m40.47s with la statis- 
tical error radius of 0'.'04. We then determined the la radius 
error circle (0'.'15) of the Chandra position of Gl by comput- 
ing the quadratic sum of the positional uncertainty for the 
X-ray source (0'.'04), the registration error between WFPC2 
and CFHT images (0'.'05), and the residuals between Chan- 
dra and CFHT alignment (0'.'13). Hence, the 95 per cent 
radius (2-dimensional) error circle will be 0'.'36. Figure 1 
shows the WFPC2 images of Gl and the 95 per cent radius 
X-ray error circle. 



4 DISCUSSION 

By utilizing Chandra, HST/WFPC2, and CFHT/MegaCam 
data, we determined the precise position of the X-ray emis- 
sion of Gl. The X-ray source, previously seen by XMM- 
Newton (see Pooley & Rappaport 2006; Kong 2007), is very 
close (~ O'.'ll) to the cluster centre. Based on the calcula- 
tion by Pooley and Rappaport (2006) , if the X-ray emission 
is from Bondi accretion of ionized cluster gas by a central 
IMBH, the X-rays should come from the central 50 milli- 
arcsecond of the cluster. However, given the uncertainty of 
the error circle, we cannot rule out that the X-ray source 
is slightly offset from the cluster core (see Figure 1). If the 
X-ray emission is not from a central IMBH, it could come 
from a luminous LMXB located within or very close to the 
core radius (0'.'21 in Ma et al. 2007; 0'.'19 in Barmby et al. 
2007) of Gl. As Kong (2007) points out, nearly half of the 
LMXBs in Galactic globular clusters are found within the 
core radius; it would not be surprising to find a luminous 
LMXB within the core of Gl. 
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Table 1. Comparison of Gl (assuming X-ray emission from Gl is from the putative IMBH), supermassive black holes, and stellar-mass 
black holes 



Source 


Mbh(Mq) 


log Lx 


log (Lx/LEdd) 


log {Ls/LEdd) 


log {Lx/Lb) 


log Ln/{Lx/LEdd) 


References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




Gl 


2 X 10*? 


36.3 


-6.11 


-4.11 


-2.00 


38.6 


This paper 


Sgr A* 


2.6 X 10^ 


33.3 - 35 


-11.22 --9.52 


-3.85 


-7.26 - -5.57 


42 - 43.7 


1, 2 


M31* 


1.4 X 10* 


35.8 - 37.3 


-10.52 - -9.00 


-4.36 


-5.45 - -3.93 


41.3 - 42.8 


3 


M32* 


2.5 X 10® 


36 


-8.52 


-5.93 


-2.83 


< 41.8 


4 


V404 Cyg 


10 


33 - 33.9 


-6.11 --5.22 






33.8 - 34.6 


5, 6 


A0620-00 


10 


30.5 


-8.62 






35.5 


7,8 


GS2000-f25 


7 


30.4 


-8.44 
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Notes. Column 1; name of the object; Column 2; Black hole mass; Column 3: X-ray luminosity (ergs~^); Column 4: X-ray to 
the Eddington luminosity ratio; Column 5: Bondi to Eddington luminosity ratio; Column 6: X-ray to the Bondi luminosity ratio; 

Column 7: ratio between radio luminosity and Lx I ^Edd- 

References. (1) Baganoff ct al. 2001; (2) Merloni et al. 2003; (3) Garcia et al. 2010; (4) Ho et al. 2003; (5) Bradley et al. 2007; 
(6) Gallo et al. 2005; (7) Kong et al. 2002b; (8) Gallo et al. 2006; (9) Garcia et al. 2001 



We estimated the probability of a l.SM© object being 
found in the X-ray error circle. The mass of I.SMq is the av- 
erage quiescent LMXB mass, estimated from the radial dis- 
tribution of 20 quiescent LMXBs in seven globular clusters 
(Heinkc ct al. 2003). Wc adopted the single- mass King model 
fit from Ma et al. (2007), which has Tc = 0.21" and also con- 
sidered the core radius value of rc = 0.19" from Barmby et 
al. (2007). We assumed that the King model profile describes 
the distribution of turnoff-mass objects; we took the turnoff 
mass to be 0.9Mq. This gives a mass ratio of g = 1.67 be- 
tween the LMXBs and the turnoff-mass stars. It is straight- 
forward to integrate the radial density profile for the LMXBs 
analytically. As a check, half of the probability is within the 
core radius for q = 1.67. We then integrated the probability 
over angle numerically, to compute the probability within 
the off-centre X-ray error circle. We found that a 1.5Mq 
LMXB would have a probability of 0.756 (rc = 0.19") and 
0.718 (vc — 0.21") of being found within our 95 per cent con- 
fidence Chandra error circle. Therefore the possibility that 
the X-ray source represents a single LMXB is significant. 

If the X-ray source is from an accreting IMBH, the ob- 
served luminosity (2 x 10'^*' crgs^^) and the spectrum would 
be consistent with a black hole in the hard state (Remillard 
& McClintock 2006). There are also two possibilities if Gl 
has an accreting IMBH. The X-ray emission could be due 
to accretion cither from cluster gas (Pooley & Rappaport 
2006), or from a companion star. 

We list in Table 1 a few examples of nearby super- 
massive black holes and dynamically confirmed stellar-mass 
black holes for comparison. All three quiescent supermas- 
sive black holes in the local universe have very low Edding- 
ton ratios of 10~^^ — 10~^. On the other hand, quiescent 
stellar-mass black holes tend to have a higher Eddington 
ratio. V404 Cyg is the most X-ray luminous stellar-mass 
black hole while A0620-00 and GS2000+25 arc the faintest 
ones (with detection). The Eddington ratio of Gl (assuming 
a mass of 2 x 10'' M©) is at least two orders of magnitude 
higher than that of quiescent supermassive black holes, but 
is in the range of stellar-mass black holes. Wc also compare 
the X-ray luminosity in units of Bondi luminosity. This is a 
better indicator of accretion luminosity than the Eddington 
ratio, as it relates the X-ray luminosity to the available mass 



transfer rate. For Gl, it is about 0.01 (see Pooley & Rap- 
paport 2006), which is substantially higher than quiescent 
supermassive black holes (see Table 1). It is an indication 
that the accretion efficiency of Gl must be high if cluster 
gas is the source of the accretion. 

In the framework of Bondi accretion, we further derive 
the Eddington-scaled Bondi accretion rate (see Table 1) and 
compare with other low luminosity supermassive black holes 
as shown in Figure 14 of Soria ct al. (2006). According to 
the standard advection-dominated accretion flow (ADAF) 
model, the X-ray luminosity of Gl should be two orders of 
magnitude lower. This implies that the putative IMBH of 
Gl is extremely efficient in accreting cluster gas. It is worth 
noting that apart from 47 Tuc (Freire et al. 2001), we do 
not have evidence of cluster gas in Galactic globular clus- 
ters (van Loon et al. 2006). Furthermore, if ionised cluster 
gas (as described in Pooley and Rappaport 2006) is the only 
source of inflow, such a high accretion rate would be diffi- 
cult to reconcile with the ADAF prediction. It is therefore 
indicative that if the X-ray emission of Gl is from an IMBH, 
it is more likely due to accretion from a companion star. 
However, based on the current Chandra data, we cannot 
distinguish between the two possible mechanisms (IMBH or 
ordinary LMXB) for generating the X-ray emission of Gl. 

We can rule out that the X-ray emission is from combi- 
nation of low luminosity sources. The total of all low lumi- 
nosity sources in 47 Tuc is only 3.7 x 10^^ ergs"'^, implying 
that the encounter rate of Gl must be > 500 times of 47 Tuc. 
However, Pooley & Rappaport (2006) estimated that the 
encounter rate of Gl is only 17 times of 47 Tuc. Hence the 
X-ray emission must come from one or a few bright sources. 
Supernova remnants may also be responsible for the X-ray 
emission but the probability of a core collapse supernova in 
a globular cluster is extremely small, as is the probability of 
a Type la supernova (e.g. Pfahl et al. 2009). Indeed, if the X- 
ray source is associated with the radio source detected with 
the VLA (Ulvestad et al. 2007), both the LMXB and super- 
nova remnant scenarios are unlikely based on the radio/X- 
ray flux ratio (sec below). The probability of positional co- 
incidence between Gl and an unrelated X-ray source is very 
small. We estimated the background AGN/foreground star 
rate within the half-mass radius of Gl (~ 1"; Barmby et al. 
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2007) using the Chandra Deep Field (Brandt et al. 2001). At 
an 0.5-2 keV flux of 1.2 x 10"^'' erg s"'^ cm~^, we expect to 
find on average 4 x 10~^ background or foreground objects 
within the half mass radius. 

Using the VLA, Ulvestad ct al. (2007) detected a radio 
source coincident with Gl at 8.4 GHz with a flux density 
of 28 ± However, the radio position is ofTset from the 

X-ray and optical centre by about I'.'S (see Figure 1) using 
our refined Chandra and HST observations. Because of the 
short baseline (3.5 km) of the VLA observations, the rms 
error of the VLA position is about 0'.'6 in each dimension. 
Therefore, the cluster core is still within the 95 per cent error 
radius (sec Figure 1). In order to confirm if the radio source 
is associated with the X-ray source, wc will require high- 
resolution and high-sensitivity radio observations using the 
Expanded VLA (EVLA) in its long baseline configuration, 
e-MERLIN, or VLBI. 

If we can confirm that the radio source is related to the 
X-ray source, this will provide a strong support that there is 
an IMBH near the centre of Gl because the radio/X-ray fiux 
ratio of Gl (~ 5 x 10~^) is several hundreds times higher 
than that of a LMXB, but is consistent with a 2 x IO^Mq 
IMBH using the relationship in the "fundamental plane" 
linking radio/X-ray fiux ratios and black hole masses (Mer- 
loni et al. 2003; see also Table 1). The radio/X-ray flux ra^- 
tio of Gl is also substantially lower than that of supernova 
remnants (~ 10~^) and low-luminosity active galactic nu- 
clei (see Ulvestad ct al. 2007 and Table 1). The radio/X-ray 
flux ratio is consistent with a pulsar wind nebula, though we 
have no evidence for pulsars with such high rotational energy 
losses having been born in old populations such as Gl. A ra^ 
dio spectrum (with EVLA or e-MERLIN) or milliarcsecond- 
resolution VLBI observations, as suggested by Ulvestad et 
al. (2007), could rule this possibility out, as could detection 
of strong X-ray variability. It is also possible that the radio 
source is not coincident with the X-ray source. In this case, 
the radio source may be a background source, a radio-loud 
X-ray binary in Gl, or a jet from the central IMBH. 
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